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1. Introduction 

Recent studies of two-dimensional Euclidean gravity have yielded a host of surprises. 
These surprises have come from lattice definitions [0-§f, continuum path-integral defini- 
tions 0-0, and topological field theory definitions [|T0HT2| of what is presumably the 
same theory 

Among the surprises from the lattice/matrix- mo del approach has been the discovery 
of theories which look very similar to the solved models of two-dimensional gravity coupled 
to matter, but for which no continuum or topological interpretation has been found. For 
example, a particularly natural family of matrix models can be derived from the integral 
over a unitary matrix. These have a double scaling limit whose exact solution is strikingly 
similar to that of the lattice gravity models, with the KdV hierarchy being replaced by the 
mKdV hierarchy (|il^l4[l . One is naturally led to ask whether these theories have equally 



simple world-sheet interpretations. 

The simplest way to attack this problem would be to find an interpretation of the 
perturbation expansion for the theory as a lattice gravity theory, possibly with matter and 
appropriate weights in the sum over graphs, and directly take the world-sheet continuum 
limit. This can be difficult since the same continuum theory can be obtained from many 
different matrix models. The important point is not the type of matrix but rather the 
structure of the saddle-point eigenvalue distribution near its points of non-analyticity. 
The mKdV scaling limits come from tuning a coupling to a critical point where two end- 
points of the eigenvalue distribution collide. This is generic for a unitary matrix, but 
is easy to arrange for a Hermitian matrix as well, for example by taking a double-well 
potential |15|,|16|]. Thus there are many possible graphical formalisms which one might 



try to interpret. Perturbation theory for the unitary matrix integral can be done by 
changing variables U = expiH; the resulting graphs have vertices of arbitrary order as 
well as fermionic lines coming from a fermionic representation of the Jacobian ||17|] . We 



will discuss this expansion further in the conclusions. Alternatively one could sum over 
the expansions about the various minima of the double- well problem. 

It is not obvious how to interpret any of these lattice theories, because we do not 
know how to take the continuum limit directly on the world-sheet. The simplest lattice 



gravity |T|-|6| is exceptional since it has the correct degrees of freedom and an everywhere 
positive measure. Nevertheless, already with non-positive measures there are subtleties in 
the continuum limit which can reproduce gravity coupled to matter |1^,|19[ . The measures 



in the present problem are too complicated for heuristic arguments. All this underscores 
the importance of developing a more precise understanding of this continuum limit. 

Since we cannot deduce the continuum theory directly from the lattice we are forced to 
compare the structure of the exact solution of the lattice theory with the structure of known 
continuum theories. With this motivation, we describe below the exact solution of the 
lattice theory in some detail. In section two we derive from the lattice the complete string 
and flow equations for two-cut models, including all even and odd perturbations in the 
potential. In section three we show how these equations can be rewritten as Virasoro-type 
constraints. This reformulation suggests the existence of a topological phase of the theory, 
discussed in section four. In section five we make some brief remarks about macroscopic 
loops in these theories. In the conclusions we attempt to draw some inferences about what 
the continuum formulation of these models might be. 

2. Multicritical Matrix Model Potentials and Physical Observables 



In this section we derive the string equation and flow formalism, generalizing 13 



14|1 , [|16f and ||2Q| . The analysis is rather technical but the results are summarized in the 
string equation (2.23), recursions (2.8), observables (2.19) and flows (2.11). 

To derive the full continuum loop equations for unitary matrix models, we will work 
with the two-cut Hermitian matrix model and study the most general scaling potential. 
The general even multicritical potential resulting in a two-cut eigenvalue distribution was 



obtained in [T(| 



V^(X)=k(m)X^+ 1 (l A2 



m = l,2,..., (2.1) 

+ 



where the subscript + means we keep the polynomial part in an expansion about infinity 
and the constant k(m) = 2 2m+1 (m + l)!(m — l)!/(2m — 1)!. 

The universality class of a potential is determined by the behaviour of the saddle-point 
eigenvalue density p(A) around the critical point. The eigenvalue density corresponding 



to V m vanishes at the critical point A = as A 2m . To obtain the most general scaling 
behaviour, we look for the potentials resulting in p(X) ~ A 2m+1 for A — > 0. Such p(A) 
are not positive definite and we will introduce them through perturbations around some 
Vm(X), in the spirit of [^TJ. The required behaviour is 

p n {\) = for |A| > 1, 



This is satisfied by 



p n (X) oc A 71 " 1 for A -> 0, 

oc(l-A) 1 / 2 for A->1~, 

oc(l + A) 1/2 for A^-l + , 
l 

dXp n (X) = 0. 
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Pn(A) = ^[A™(l-A 2 ) 3 / 2 ], n=l,2,... 



(2.2) 



and the corresponding potential is 

V n = v(n)N c ^X n (X 2 -lf /2 \ + , 

with v(n) and c(n) constants to be fixed. 

The scaling ansatz appropriate to the potential V m is 

R n = r c + (-ira 1 / m f(x) + --- 1 

S n = s c + (-l) n a 1 / m ig(x) + ---, (2.3) 
^ = l-a 2 (x-p), Na 2+1 / m = l. 

s c is of the order of odd perturbations, which are supressed by a power of N . The or- 
thonormal polynomials V n (X) tend to two scaling functions, depending on parity: 

(-l) n P 2 n(a 1/m A) = a 1/2m p + (x. A), 

(2.4) 

(-l) n V 2n +i(a 1/m X) = a^ 2m p-(x-a^ m ,X). 

From (2^) and ( |2.4f) it follows that the multiplication by the eigenvalue A = a 1 / 771 A is 
represented on \1/ = ^+ J by a 2 x 2 matrix (after a rescaling of / and g and a rotation 
of*): 

\J 3 V = (D + q)V, (2.5) 
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where D = d/dx, q = f(x)Ji + g(x)ij%, and Ji are matrices satisfying [J7i, Jj] = iSijkJk- 
Finally, redefining — > we obtain (dropping tilde from A) 



We will need the resolvent 



R 



1 



(D + q)J a -$ 



A • 



In their work on the resolvents of matrix differential operators p2|| Gelfand and Dikii show 
that R satisfies 

J 3 R f = [J 3 R,q-XJ 3 }, (2.6) 
and that it has an asymptotic expansion 



R(x,X) = J2 R kX~ k - 

k=0 

From ( |2.6| ) it is clear that, up to constants, the most general expansion of R is 

oo 

J 3 ^ = ^grad^ fc _ 1 A- fc , 



(2.7) 



k=0 



where grad Hk = —JiGk—ifoFk + JzHk- Plugging (|2.7| ) into (|2.6|) results in the recursion 
relations determining Gk, Fj~, Hj-. Setting G_i = F_ x = 0, = 1, we obtain 

Fk+i = G' k + gH k 

G k+1 = F' k + fH k (2.8) 
H'k = dGk - fFk, 

for all k > 0. The first few are 

-Fb = 9-, Go = f, 



Fi = /', 



#n = 0, 



^i = ^(^ 2 -/ ; 



(2.9) 



r)> H 2 =gf'-fg', 



*3 = /"' + |/V - f 2 ), G 3 = g"> + ~-g\g* - f 

Hs = gg" - \{g') 2 - ff" + \{f) 2 + 3 -(g 2 - f 2 ) 2 . 



This determines the expansion of the resolvent 

oo 
k=Q 

There is an SO(l, 1) symmetry in our system coming from possible choices in the original 
definition of q. Observables will be invariant under this symmetry. 

A related derivation of (|2.8| ) would start from two Hamiltonian structures defined in 
||23|| for the system with the Lax operator 

L = D + q-XJ 3 . (2.10) 

For the flow in the coupling t k: with dimension k, their compatibility gives 

J-g = [grad H k , J z ) = [grad # fc _ l5 D + q], (2.11) 
atk 

which is easily seen to result in (|2.8|) . 

We are now ready to start determining the observables. We will calculate the 2-point 
functions of the puncture operator P with the scaling operators a n corresponding to the 



perturbations by V n . We will follow the presentation of [21], beginning with the one-point 
function 



dF 



dt 



it 



v(n)N^ I ^-Y{V 2 - lf /2 \ + tr( -^IlA (2.12) 

Jc 2m \y - A J 



where 

m 1 ifO<i,i<JV-l, 

{ Nhj \ otherwise, 
and C surrounds all of the eigenvalues of A. By deforming C to a circle at infinity one can 
see that the subscript + in ( |2.12j ) can be dropped. To obtain a two-point function, we 
study a variation 5dF/dt n and use 

Str (-^ U N) = ltr(^—[S(6V(A)),U N ] 



y-A J 2 \y-A [ 
S(B) ij = e(i-j)B ij , 

= sign of j, j ^ 0, 

= 0, j = 0. 
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We choose the simplest even operator as the puncture operator. The corresponding vari- 
ation of the potential is SV(A) = N Cp A 2 , and we obtain (to leading order in N): 



= v(n)r c N Cp+c ^ j ^^(y 2 ~ !) 3/ 



1 \ I 1 

+ 



(2.13) 



The scaled quantities are 



y 


= a 1 / m A, 


A 


= a 1 / m 4(D + q)J 3 , 


\2n) 


= (-i) n a 1/2m \x,+), 


\2n + l) 


= (-l) n a 1/2m \x-a 1/m ,-), 



(2.14) 



and ( |2.13| ) becomes (again to leading order) 
d 2 F 



dxdt n ^) r i^ +1),m ^ P+c{n) I +\R\*, +) + <*, -1*1*. -> W 



, C Z7TZ - (2.15) 

= v(n)^-a {n+1)/rn N Cp+c{n) H n (^). 

The simplest non-trivial case in this formalism is m = 1 potential ( ^.1| ) with simplest 
odd perturbation []20[ 

V = Vi + Naetr(<f>). (2.16) 

(Note that, in the context of the unitary matrix model as 2D lattice Yang-Mills theory this 
coupling corresponds to a discretized version of logdett/ mod27r, and hence corresponds to 
a lattice version of the theta angle ||24|1 .) Using ( p.3|) with s c = — ae/16 in the recursion 
relations gives string equations^ 

= x f + eg' + 2(f" + }-f(g 2 -f 2 )), 

7 (2-17) 
= X g + ef' + 2(g" + -g(g 2 -f 2 )). 



1 To obtain these string equations we rescale e, / and g and shift x. In order to compare with 
|2.9| ) the rescalings of / and g have to be the same as the ones used to obtain ( [2.5D . 



The specific heat is 

(pp> = -i + s + i (92 - /2) - (2 - 18) 

The first two terms are non-universal and the universal piece is given by Hi/2. Therefore, 
we can identify x = t±, o\ = P. This identification shows that, in order to obtain a finite 
result in (|2.15| ), we need c(n) = n/(2m + 1). II Finally, for later convenience, we choose 
v(n) = l/r c for every n. In conclusion, in the scaling limit of the two-cut matrix model 
we have obtained the following observables 

d 2 F 



dxdt 



(Pa n ) = \H n . (2.19) 



n 



The odd flows F%k and G2k have terms with no derivatives. However, an insertion of 
an odd operator is still suppressed by 1/N. 
To derive the string equations we will use 

n = RT [ V'VnVn^ = I ^V'(y) (-*-) , (2.20) 

J Jc 2m \y- A Jn,n-l 

2m 

V = Nk(m)y 2m (y 2 - l) 1 /^ + ^ tk L N ^ y h-i {y 2 _ 1)3 /2 |+ + . . . j (2 21) 

In (|2.21| ) + ■ • ■ denotes terms of higher order in a upon introduction of scaling variables 
through n = N(l — a 2 x) and (|2.14|) . 

We want to replace the y-integral along C, a contour surrounding all of the eigenvalues 
of A, by an integral over A, i.e. over infinitesimally small y's. An obstacle to passing to a A- 
integral comes from the A-eigenvalues outside of the critical region. Their contributions will 
not have scaling dependence on coupling constants, but can be large (O(N)) nonetheless. 

In order to be able to neglect consistently eigenvalues of A at large y, instead of scaling 
( p.20|) , we will scale its derivative with respect to ti, I ^ 1. Namely, the ti dependence of 
V appears only for y very small (C(a 1 / m )), while for y ~ 0(1) it is suppressed by a power 
of (large) N. The result is 

°=IT { ^T,kt k X k -\T)(x,±\f](-iJ2G r . 1 -J 1 F r _ 1 +H r _ 1 /4)X- r \x 1 T), (2.22) 



k r=0 



2 The factor of Na in front of tr((f)) in (|2.16|) is precisely iV 2 / 3 appearing in V2. We use the 
fact that the scaling limit of a generic odd (even) potential like tr(4>) (tr(0 2 )) is given by V2 (Vi). 



where the upper (lower) sign corresponds to n even (odd). After doing the integral and 
taking the 1, 2 (2, 1) element of the resulting matrix, we have 



= ^kt k G k -i = Sg, 

k>l 

= J2 kt kF k -i=S f . 



(2.23) 



fc>i 



The integration constant, which is independent of all ti, I ^ 1, can be seen to be zero by 
comparison with (|2.17| ). (It was shown in |E3] that by modifying the action to include 
boundaries, or "quark" degrees of freedom one can introduce a nonzero constant on the 
left hand side of (ggg).) 

The string equation generalizes ( 2.17 ). Moreover, it may be equivalently written as a 
flatness condition, generalizing that in ||16|| : 



d 

L 'Tx- M 







(2.24) 



where 



k>i 



(2.25) 



M k = grad H k _ x + A grad H k _ 2 + ■■■ + A fe_1 grad H + X K J 3l 
and Mq = J3. M is chosen to make the left hand side of ( |2.24|) independent of A, as 



described in 23 



The flatness condition (|2.24j ) can be interpreted as the isomonodromic deformation 
condition for the monodromy defined by a solution of 

d 



dX 



— M vp = 0. 



(2.26) 



Furthermore, the partition function is again the isomonodromic r-function. The proof of 
these statements follows the steps already outlined in |26|,|16|]. i Moreover, the solutions 



The form of recursion relations useful in the proof is ( |2.8| ) with the third equation replaced 



by 



1 

Hk + l = — ^ ^ (F r Fk-r — G r Gk-r — H r Hk-r)- 



r=0 
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to the string equations chosen by the matrix model will in general lead to nontrivial Stokes 



data. As described in [p6[ this implies that the r function does not correspond to a point 
in the Segal- Wilson Grassmannian, although it does correspond to a point in the Sato 
Grassmannian. 

3. Loop Equations. 

The solutions of one-matrix models satisfy Schwinger-Dyson equations which can be 
derived by variation of the functional integral ^7|. These have a double scaling limit 



|2£|^9| which in the KdV-type systems can be rewritten as Virasoro constraints on the 



partition function |29] [|30 



Since the same Schwinger-Dyson equations apply (before taking the continuum limit) 
to the two-cut Hermitian matrix model, and since the unitary matrix model has analogous 
Schwinger-Dyson equations, one would expect continuum loop equations to exist in this 
case as well. 

Previously we had derived a subset of these loop equations (with two non-zero cou- 



plings) by directly scaling the matrix model Schwinger-Dyson equations [pl[ . Here we 
derive the complete loop equations from the string equation/mKdV formalism. The argu- 
ment is similar to that used in |3(J for the KdV theories. Using the flow equations, the 



string equation can be reformulated as a differential constraint on the partition function 
as a function of the couplings. In principle we can then derive a series of constraints as 
follows: the same recursion relation which takes the n'th flow to the n + l'st can be applied 
to the n'th constraint to find a new constraint. 

Since the constraint operators all annihilate the partition function, they form a left 
ideal. Therefore commutators of constraints are also constraints. Guided by our expecta- 
tion that the constraints form a Virasoro algebra, we will find a minimal set of constraints 
which generate this algebra. The completeness of the resulting algebra will then follow if 
we can argue that there is a unique solution of the resulting constraints. This will be true 
in a certain sense to be discussed below. 

The scaling operator corresponding to the coupling t n will be called a n . a\ will also 
be called P, though one should probably think of the P as standing for "primary" and not 

10 



"puncture" for reasons explained below. The string equation is now Sf = S g = 0, and our 
constraints will be derived from combinations of these. 
We start with 

= gS g -fS f 
= Y,kt k {a k ^PP) (3i) 

= D 2 ^2ktk(a k -i), 

k>2 

which we interpret as the second derivative of an L_i constraint. The explicit dependence 
on t\ has cancelled out. 

The recursion relations imply that gG k +i — fF k+1 = gF k — fG' k , so we can get another 
equation from 

= gS' f - fS' g 

= J2kt k (a k PP) + (g 2 -f 2 ) 
k 

= J2ktk(*kPP)+2(PP) (3 ' 2) 

k 

= D 2 ^kt k {a k ) 

k 

the second derivative of an L constraint. 

We can continue in this vein using the recursions to construct t k d/dt k + n in terms 
of differential polynomials of the string equations. We will skip L\ (which is similar) 
and proceed to the Li constraint, since it will generate the complete algebra. Using the 
recursions three times, we have 

k 

= \ Yl kt k(gG k+2 - fF k+2 ) 
1 k 

= 2 kt k(gF k +i - fG'k+i) (3.3) 

k 

= \Y, kt «(9G k ' - fF>> + (g 2 - f)H' k + (gg' - ff)H k ). 

k 

= IJ2 kt k{9F'k-i - fG' k n _ x + X fc _! + (g 2 - f 2 )H' k + (gg' - ff)H k ). 



2 

k 
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with = g(fHk-i)"—f(gHk-i)"- The X piece is identically zero by the L_i constraint 

( [3J]) and its derivatives: 

D n ktkWk-i) = ktkD^H^. (3.4) 

k>2 k 



Using S'f = J2 ktkF^ + 3g" = (resp. for g), this i 



is 



h = -%g" - ff") + 2j2kt k (2{PP)(a k PP) + (PPP)(a k P)) (3.5) 



k 

and using derivatives of the L Q constraint 



h = -8(PP} 2 - 2(PPP)(P) - \{gg" - ff"). (3.6) 
The other terms we expect in an L 2 constraint are 

(PPPP) = \ (gg" - ff" + {g'f - (f) 2 ) , (3-7) 

foP) = \W - ff" - \{ 9 '? + 1(f) 2 + l(9 2 ~ f 2 ) 2 ). (3-8) 
These combine nicely to give 

= J2 kt k (a k+2 PP) + 2(a 3 P) + (PPPP) + 2(PP} 2 + 2(PPP)(P). (3.9) 

k 

Defining Z = e F , this is the second derivative of the expected constraint L 2 Z = with 

i *=E"»«ib + M- (3 ' 10) 

k>l K+z 1 

Note that to get to ( |3.5| ) we assumed that the undifferentiated L_i and Lq constraints 
were satisfied, i.e. that they contained no integration constant. Under this assumption, 
we can now write the result as Virasoro constraints for a two-dimensional free massless 
boson. Define 

d<p( Z ) = kzk ~ H * + E z ~ x ~ k w> (3 - n) 
k>i k>i k 

an untwisted boson; then the constraints are the modes L n for n > — 1 of the usual 
untwisted stress-tensor (dtp) 2 with zero ground-state energy. 

12 



The integration constants in the constraints must be compatible with the Virasoro 
algebra, otherwise the commutator of two constraints would generate new constraints 
which would force Z to be trivial. An obvious possibility would be to include the zero 
mode of the boson as well, 



q would be a new parameter of the theory, not associated with any operators or flows. In 
the "usual" phase of the theory with all odd couplings set to zero, it is clear that q = 0, 
since all the other terms in the L_i constraint are zero, but we will find a use for non-zero 
values shortly. 

It should also be possible to derive these loop equations directly by double scaling 
the finite iV Schwinger-Dyson equations. In particular, the theory of the boson ( |3.12j ) 
should emerge from the formulation of the lattice theory using second quantized fermions. 
Unfortunately, the currently available derivations of this "field theory on the spectral 
curve," are merely heuristic. 

4. A Topological Phase 

In the series of multicritical models including pure gravity, described by the string 
equation 



the simplest point is not pure gravity but rather the "pure topological" theory with t\ = 1 
and t n = for n > 1. This theory has a conserved charge ("ghost number") with a distinct 
background charge for each genus surface, so a specific correlation function can be non-zero 
for at most one genus. In particular ((rococo) = w' = 1 is non-zero only on the sphere. 

With purely even couplings, there is no analogous model in the mKdV hierarchy, but 
with odd couplings there is a precisely analogous model with t2 — 1 and t n = for n > 2. 
The string equation is then 




(3.12) 



X = t\U + t2(u" + 3w 2 ) + . • • , 



(4.1) 



xg + f = 



(4.2) 



xf + g' = 0. 
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with solution 

9 + 1 = 0^/* 

(4-3) 

9-f = C 2 e* 2 '\ 

C\ and C2 are constants of integration. Although these solutions are trancendental, phys- 
ical quantities are not: for example, 

(°i°i) = \{9 2 -f 2 ) = \c 1 C 2 , (4.4) 

so this model appears to have "7 string = 0" . (The Baker functions, which are Airy functions 
in topological gravity, are parabolic cylinder functions in this theory.) 

The SO (1, 1) symmetry insures that all correlation functions in this model are poly- 
nomial in x. Assign the combinations f + g and / — g SO(l, 1) charge 1 and —1; then 
correlation functions must be neutral. All charge-neutral polynomials in these functions 
are constant, while appearances of d/dx will produce positive powers of x. 

The structure of the topological phase is seen most clearly in the loop equations. Just 
as for topological gravity, expanding the loop equations around these couplings gives a 
series of recursion relations, one for each operator in the theory, expressing an operator 
insertion in terms of correlators of lower total degree. For example, the L_i constraint 
gives us 



2 



(where a = 0) for all cases except one: 



(a 1 a 1 ) = -\. (4.6) 

To get non-trivial answers in the topological phase, we must take the parameter q = —C1C2 
non-zero. 

Just as for the KdV models, the existence of a phase in which all operator insertions 
satisfy recursion relations, means that the expansion of the partition function to all orders 
is uniquely determined by the constraints. They can therefore have only one analytic 
solution. Unlike the KdV models, since the parameter q is not a coupling of the model, 
we cannot flow from this topological phase to the phases with q = 0. This is probably 
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connected with the difficulty of defining this pure topological phase in the matrix model 
- our odd perturbations of the potential ( p.2|) really only make sense as perturbations of 
even potentials. On the other hand the structure of the pure topological phase seems very 
clear from the solution, so we will not let this stop us from studying it. 

The ghost number conservation law follows most directly from the loop equations. 
If we write them in terms of our boson d(p(z), they are homogeneous in z, so the ghost 
number assignment of each operator and coupling must be just the associated power of z 
in a mode expansion, up to a possible constant shift. We can use this constant shift to 
make the ghost number of one coupling of our choice zero. If we choose t2, setting t2 = 1 
will not break ghost number conservation. 

For topological gravity this is the end of the story. The special role of £2 in this theory 
is that it is the only choice of non-zero coupling for which all operator insertions have an 
associated recursion relation. In the present case, both and q must be set non-zero, and 
they have different ghost number, q having that of "to" ■ A simple example of a correlation 
function which shows that q must be assigned ghost number is 

(a 3 a 4 ) = -48AV - 144Ag 3 x + 168AVx 3 - 24AV 5 (4.7) 

where A counts derivatives (powers of 1/N). 

We can maintain a ghost number conservation law if we take q non-zero but treat it 
as the coupling of a new operator in the theory. If the operator a n is given ghost number 
2 — n, then q will couple to an operator Q with ghost number 2. Now the basic correlation 
function on the sphere is 

(wiQ) = ~\ (4-8) 

Q is an unusual operator in that it does not have descendants. The ghost number 
counting would suggest that it could be renamed 00 and that the other operators are its 
descendants, but this is misleading. It clearly appears in a very different way in the string 
equation; in the loop equations, it is different in that its conjugate d/dq does not appear. 
The necessity to include this operator is the most serious difference with the topological 



gravity coupled to matter of [10 
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5. Macroscopic Loops for Unitary-Matrix Models 

It has recently been shown that the study of macroscopic loop amplitudes |32] is 
particularly well-suited for the comparison of continuum-Liouville and matrix model for- 



mulations of 2D quantum gravity [|33| |34|] . In particular, the Wheeler-DeWitt equation 
plays a central role in making such comparisons. In this appendix we derive analogous 
macroscopic loop amplitudes for the unitary-matrix model, and show that they too obey 
a Wheeler-DeWitt-like equation. This equation gives a hint about the proper worldsheet 
interpretation of these theories. 

The macroscopic loop operator is defined to be ^e~ iL ^ in the fermionic formulation, 
where \1/ is a two-component fermi field (obtained from scaling the even and odd indexed 
orthogonal polynomials) and L is the Lax operator representing multiplication by A. The 
one-point function is therefore 



X 




w(£)}= / dx(x\tile \J y " / \\x) (5.1 



We obtain the genus zero approximation by evaluating the matrix element with the WKB 
approximation to get 



( W (e)) 




f - 9 -ip 



(5.2) 



+ e Vs 2 -/ 2 -j 

p oo rOO 

= 2 dx dp cos(V:P 2 - #i) 

J ti J — oo 

POO 

= -tt / dxh\ /2 Ji(lh\' 2 ) 
Jtl 

where hi = —H\ is positive (for large x) and J\ is a Bessel function. Using the flow 
equations, specialized to the case of genus zero, we therefore obtain the wavefunctions of 
the even flow operators: 



/OO or 
dxh^^JoiihY 2 ) 

l>00 

= -nl / dyy^Miy 1 / 2 ) 

Jhi 



(5.3) 
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The operators for the odd flows have zero wavefunction at tree level. The integral in 
( |5.3|) only converges in a distributional sense. We may imagine regularizing the original 
macroscopic loop operator so as to obtain an additional factor of e~ £ ^ in (|5.3|) . It is then 
easy to show that (|5l^) and subsequent formulae below possess a finite unambiguous limit 
as e — > + . In particular, using the identities ^-J„(z) = J v+ i(z) + J v -\{z) and 

/>oo 

lim / d y e-^{ y - if" 1 J Q (ay/y) = T(k)2 k (-a)- k J k (a) (5.4) 

we may prove (as in [[J3|) the existence of a linear, upper-triangular, analytic change of 
operators a v — > a v which have the wavefunctions: 

1> v (£) = (a u w(£)) = h\ /2 J u {lh\ /2 ) (5.5) 

for v = 2k- 1, with k = 1,2, .... 

The wavefunctions of the a operators satisfy the Bessel equation, which, in the one- 
cut phase has been interpreted as the Wheeler-DeWitt equation for gravity coupled to 
conformal matter in the minisuperspace approximation. Assuming a similar interpretation 
holds in this case we see that the Bessel equation should be interpreted as 



ip„ = (5.6) 



where 7, Q, [i are parameters in the Liouville theory, following the conventions of |3q] , and 
A is the dimension of the matter field being dressed. It follows that the cosmological 
constant must be considered to be negative. Accordingly, for a convergent path integral 
the worldsheet should have a Minkowskian signature [BBJ. Moreover, it follows from ( |5.(j| ) 
that the central charge and a subset of dimensions of the matter theory must satisfy 

8 ,1-c 

+ Aj 

(5.7) 



» 2 = v^ +A) 



7 2 = ^(13-c- v / (l-c)(25-c)) 
o 

where v is any odd integer. 

We can now use ( ^75|) ( |5.7|) to learn about the continuum theory. Let us assume that the 



structure of wavefunctions ( p~5| ) is exactly analogous to that described in |33| . In particular, 
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we assume that the multicritical phases of the matrix model defined by k = 1,2, 3, . . . are 
characterized by the identification of the wavefunction of the cosmological constant with 



ip Uo for vq = 2k — 1. Plugging into (|5.7| ) we find a series of matter theories with central 



charges 



6^ 



c = i _ (5.8) 

Again using (|5.7|) we discover that the are wavefunctions of operators whose flatspace 
dimensions are 



v- 



The above critical exponents are consistent with the identification of the matter theory 
as a special case of the 0{n) model [|37|]- In particular, the 0(n) model for n = — 2cos7rgr is 
thought |38j to have a hierarchy of multicritical points, for a fixed n, with central charges 

6(S-1) 2 



1 



9 



(5.10) 



Moreover, the spectrum of the model includes operators with weights given by |39| 



(5.11) 



for v G 7L. The operator algebra respects the TLi grading defined by the parity of v. 
It is possible to identify the TL<i eigenspaces as Ramond and Neveu-Schwarz sectors of 
a conformal field theory |[40[| . In view of these results [|38]] [[39]] we may identify g = 
uq + 1, an even integer, so that n = —2. Only the "Ramond operators" with v odd have 
nonvanishing wavefunctions at genus zero, but it is natural to identify all the a v as the 
gravitationally dressed versions of the 0(n = —2) operators with dimensions ( |5.11|) . 



6. No Conclusions 

In the previous sections we have described in detail the solution of the unitary-matrix 
model, including string equations, flows, Virasoro constraints and macroscopic loop ampli- 
tudes. Returning to the questions mentioned in the introduction, we can now make some 
comparisons with candidate continuum formulations. We will discuss three candidate the- 
ories below. Unfortunately, in all three cases the available evidence remains inconclusive. 
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The first candidate is 2D supergravity. The multicritical points would correspond 
to supergravity coupled to the (2, 4k) superminimal matter theories, as already discussed 
in 



5TJ. Evidence for this is (1) the structure of even and odd scaling operators and their 



dimensions match the (2, 4k) series of super-minimal matter; (2) correlators of even (Neveu- 
Schwarz) operators agree with those of bosonic gravity at tree level [[|l|; (3) there is no 
space-time supersymmetry or Grassmann couplings, as one might expect for supergravity 
coupled to low c minimal models with no possibility of GSO projection. 

Clearly the above is not very compelling evidence. One would like to have, for ex- 
ample, tree level Ramond correlators or an understanding of the super- WdW equation to 
make a more definite statement. Higher genus correlators could also settle the question 
but problems associated with sums over spin structures make these difficult to calculate. 
In particular, the calculation of the one-loop partition function from the continuum the- 
ory could be revealing. The contribution for the even spin structures has been carefully 



determined in |[42|| . There is a further point which is rather difficult to reconcile with this 
interpretation. The two-cut model in one dimension is the same to all orders in pertur- 



bation theory as the one-cut model [J3| ||44|| . This would be rather surprising for c = 1 
matter coupled to supergravity. Nevertheless, the supergravity interpretation has not yet 
been definitely ruled out. 

The second candidate is dense polymers or the 0{n = —2) fermionic model, (these 
two being intimately related non-unitary matter theories with c = —2 coupled to 

Lorentzian metric gravity. The principle evidence for this was described in section five 
and follows from assuming that the macroscopic loop amplitudes satisfy the Wheeler- 
DeWitt equation of 2D gravity. A further hint that this is the correct interpretation comes 
from examining in detail the perturbation series of the simplest unitary-matrix model 



using the action in []17 1 . The Feynman perturbation series can be written as a sum over 
surfaces together with a sum over self-avoiding loop configurations (the fermion loops) on 
those surfaces. Unfortunately, the Boltzman weights for the loops are similar to, but not 
precisely equal to those of the 0(n) model and it is not clear if the difference is important at 
criticality. Moreover, in this interpretation it is also rather mysterious that the partition 



functions at every genus should be positive [|13|]. Finally, the ^2-odd loops should be 
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calculated and compared with the NS sector of the polymer problem. 

The third candidate is a candidate for the topological phase discussed in section four. 
The similarity of the structure of the topological phase to the KdV-type systems strongly 
suggests that the theory has a topological interpretation. Ghost number counting rules out 
the obvious possibility that it is Osp(2\l) topological gravity. Another idea follows from 



recalling the identification of t% with the "Yang-Mills 9 angle" described below (|2.16| ). In 
the continuum this operator becomes simply J trF which could serve as an action for a 
topological Yang-Mills matter theory, which could then be coupled to topological gravity 



[fR)| | IH ] . However, there are some troubles with any candidate topological matter 
theory. If we have a theory of topological gravity coupled to topological matter, then, as 
shown by Dijkgraaf and Witten there must be a "puncture operator" P which satisfies a 
"puncture equation" expressing an insertion of P in terms of correlation functions of lower 
degree. The natural candidate for the puncture equation would be the L_i constraint. 
However, this interpretation would make all operators a n descendants of a single primary, 
o"i, which seems incompatible with the uniqueness of the known TFT with a single primary. 
In addition, ghost number conservation required us to introduce a new operator Q with 
no descendants, which does not fit at all into the topological gravity framework. 

Although it is possible that a different interpretation of the model exists, which is 
more compatible with the original topological gravity framework, we did not find such a 
description and suspect that it does not exist. If this is true, clearly it would be very 
interesting to find the correct world-sheet topological description. 
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